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ABSTRACT 

Discoveries of rotating radio transients and fast radio bursts (FRBs) in pulsar surveys 
suggest that more of such transient sources await discovery in archival data sets. 
Here we report on a single-pulse search for dispersed radio bursts over a wide range 
of Galactic latitudes (|6| < 60°) in data previously searched for periodic sources by 
Burgay et al. We re-detected 20 of the 42 pulsars reported by Burgay et al. and 
one rotating radio transient reported by Burke-Spolaor. No FRBs were discovered in 
this survey. Taking into account this result, and other recent surveys at Parkes, we 
corrected for detection sensitivities based on the search software used in the analyses 
and the different backends used in these surveys and find that the all-sky FRB event 
rate for sources with a fluence above 4.0 Jy ms at 1.4 GHz to be 7?. = 4.4^3'^ x 
10^ FRBs day“^ sky“^, where the uncertainties represent a 99% confidence interval. 
While this rate is lower than inferred from previous studies, as we demonstrate, this 
combined event rate is consistent with the results of all systematic FRB searches 
at Parkes to date and does not require the need to postulate a dearth of FRBs at 
intermediate latitudes. 
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1 INTRODUCTION 


observed in the 1-2 GHz band. To date, 11 FRBs 



have been published^ 

(iLorimer et al. 

20071: iKeane et al. 

The radio transient sky contains a number of known and 

2OI2I: Thornton et al.r 

2OI3I: Thornton 

2OI3I: Snitler et al. 

potential classes of sources which emit on timescales rang- 

2OI4I: Burke-Snolaor & Bannistei]l2014l:lPetroff et al.ll2015al: 

ing from nanoseconds to years. These classes include the 

Ravi et al.ll2014l'). Although, most of these were discovered 


high energy particles, magnetars, y-ray bursts, maser flares, 
active galactic nuclei, radio supernovae, pulsars, annihilat¬ 
ing black holes, and transmissions from extraterrestr i al civ¬ 
ilizat ions (for a review, see, e.g., Cordes et all l2004l . [Laziol 


I2OI2I I. Recently IPietka et al.l ll2015[ ) have analyzed the vari¬ 
ability timescales of radio flares from nearly 90 different ob¬ 
jects/events varying from flare stars to supermassive black 
holes in active galactic nuclei and have demonstrated that 
variability timescales could be used as an early diagnostic of 
source class in future radio transient surveys. 

Within the last decade, as part of the analyses of 
pulsar surveys, a new transient phenomenon — known as 
fast radio bursts (FRBs) — has been identified. FRBs 
are categorized as short duration (few ms) bursts that 
are non-repeating with dispersion measures higher than 
the Galactic di spersion measu r e (D M) expected along 
that direction dCordes fc Laziol l2002tl . The DMs range 
from 300-1600 cm“® pc and have so far only been 


FRB 131104 (iRavi et al.ll2014l l are real-time discoveries with 
a transient pipeline developed at the Parkes telescope. De¬ 
spite a significant amount of follow-up observatio ns which 
closely followed FRB 140514 (iPetroff L aklf^lSal l. and in¬ 
spection of archival surveys for other FRBs, none of the 
FRBs currently known are associated with any source coun¬ 
terpart at other wavelengths. 

As remarked by many authors, the high DMs of 
FRBs are suggestive of an extragalactic origin. Among 
the proposed extragalactic sou rce populations for FRBs 
are annihilating black holes (|Keane et al.l I2OI2I I. flar- 


(IPodov & Postnov 

I2OI3I'). binary white 

llKashivama et al. 

2013|b binarv neutron 


star mergers (iTotanil l2013l i. collapsing neutron stars 


^ For an updated list, see http://astro.phys.wvu.edu/FRBs 
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dFalcke fc Rezzoi5l2014 ), and neut ron star-black hole merg¬ 
ers i Llpunov fc Pruzhinskav3l20l4l . 

In view of the lack of any extragalactic counterparts 
identified s o far, a number of other scenarios remain equally 
intriguing. iLoeb et al.l (l2014h considered the case of nearby 
Galactic flaring main-sequence stars within 1 kpc as the 
sources of FRBs. They propose that the excess dispersion 
com es from propagation thro ugh the stellar corona. How¬ 
ever [Lua^^^^Goldrei^ ll2014ll argued that the free-free ab¬ 
sorption would conceal any radio signal emitted from below 
the corona. Also, FRBs exhibit quadratic dispersion curves 
that are consistent with the assumption of weak dispersion 
in a low density plasma. So, for the above model, in which 
the dispersion is concentrated in a relatively high density re¬ 
gion, the qnadratic dispersion approximation breaks down as 
the plasma frequency is comparable to the propa gation fre¬ 
quen cy, posing s ignificant pr o blems for this model JPennisorJ 
|20141 . However. iMaoz et al.l ([201^ have recently fonnd pos¬ 
sible flare stars in additional FRB fields using time-domain 
optical photometry and spectroscopy. The authors have eval¬ 
uated the chance probabilities of these possible associations 
to be in the range 0.1% to 9%. FRB 140514 was discovered 
in the radio follow-up observations of FRB 110220 , three 
years apart within the same radio beam. iMaoz et al.l (l2015l l 
also claim that these two FRBs are from the same repeating 
source with 99% confidence and are consistent with the flare- 
star scenario with a varying plasma blanket between bursts. 
More FRB detections in general are necessary to confirm or 
refute a Galactic origin of FRBs. 

Although FRBs have so far been observed over a range 
of Galactic latitudes, their true distribution on the sky re¬ 
mains unclear. Recently, based on an analysis of Parkes 
High Time Resolution Univers e (HTRU) mid-latitude sur¬ 
vey data, IPetroff et al.l ll2014h proposed that there is a 
deficit of FRBs at intermediate latitudes. Further analy¬ 
ses of archival and current s urveys are required to inves¬ 
tigate this issue. In addition, iKeane fc PetrofS (l2015ll have 
assessed the commonly used search algorithms used for 
FRB searches which impact the FRB sensitivities in indi¬ 
vidual surveys. Motivated by these results, in this paper 
we pres ent a single-pulse search which is sensitive to both 
RRATs llMcLaughlin et ^l2006ll and FRBs on archival data 
previously searched for pulsars at high Galactic latitudes by 
iBurgav et al.l ll2006ll . In § 2 we outline the basic methodology 
behind the search and the detection criteria. Our results and 
single-pulse statistics of known pulsars re-detected in this 
search are discussed in § 3 . In § 4 we use our non-detection 
of FRBs to constrain their all-sky event rate. Finally, in § 5 
we summarize our results and present our conclusions. 


2 SEARCH METHODS AND ANALYSIS 
2.1 Survey parameters 

The Parkes high-latitude (PH) survey dBurgav et akllioo^ l 
was designed to find millisecond pulsars and exotic bina¬ 
ries which migrate away from the Galactic plane. A total 
of 42 pulsars were detected in this survey, of which seven 
were millisecond pulsa rs and 18 wer e new discoveries. The 
analysis of the data bv lBurgav et al.l (l2006tl was carried out 
using the standard periodicity search methods to find pe¬ 
riodic signals. However no single-pulse searches have been 


published on these data. The single-pulse search method is 
very effective in detecting sporadic sources like some pul¬ 
sars (e.g. nulling pulsars) and RRATs in the time domain 
since these might not be detectable in the standard periodic¬ 
ity searches (ICordes fc McLaughlin|[2003l l. Moreover, FRBs 
can of course only be detected through single-pulse searches 
(see below § 2.2.2). 

The PH survey covered a strip of the sky enclosed by 
Galactic longitudes 220 I ^ 260 and Galactic latitudes 
|6| 60 corresponding to a total sky coverage of 3588 deg^ 

in 475 hours. The survey began in November 2000 and 
ended in December 2003 and made use of the 13-beam re¬ 
ceiver on the Parkes 64-m radio telescope. Data were col¬ 
lected simultaneously by 13 beams at a central frequency of 
1374 MHz with 96 frequency channels, each 3 MHz wide. 
Each of the 6456 pointings was observed for 265 s. For more 
details about the re ceiver system and data acquisition, see 
iBurgav et al.l (l2006tl . 

2.2 Single-pulse search method 

Each of these beams from the survey was processed inde¬ 
pendently using the SIGPROC software packagfl The steps 
included in our analysis are: 

(i) dedisperse the raw data file at a range of trial DM values 
and remove radio frequency interference (RFI) at zero DM; 

(ii) search for individual pulses in the time series above 
signal-to-noise (S/N) of five and with different widths; 

(iii) apply the detection criteria to filter in terms of DM, 
S/N, and number of beams; 

(iv) manually inspect the resulting diagnostic plots. 

We describe each step in detail below, and give the appro¬ 
priate SIGPROC modules used. 


2.2.1 Dedispersion 


Radio signals are affected by interstellar dispersion, and 
as a result, the higher frequencies of the signal traveling 
faster through the interstellar medium arrive earlier than 
their lower frequency counterparts. The time delay between 
the t wo frequencies fi and / 2 , (see, e.g.. iLorimer fc Krameil 
l2004ll is 


At ~ 4150s 


/i 

MHz 


-2 



/ DM \ 
\cm~3 pc J 


( 1 ) 


where DM is the integrated number density of free electrons 
along the line of sight. This dispersion allows us to distin¬ 
guish between astrophysical and terrestrial signals. It causes 
a quadratic sweep across the band and may be removed by 
appropriately shifting the frequency channels. Each time se¬ 
ries was dedispersed over a range of trial DM values using the 
dedisperse routine in SIGPROC. For this analysis, we have 
searched DMs in a range from 0 — 10^ cm“® pc. The wide 
range of DM makes the search sensitive to events that are 
highly dispersed. The trial DM step sizes used in the anal¬ 
ysis were calculated u sing an algo rithm, dedisperse_all, 
originally described by iLevinI (120121 ') which accounts for the 


^ http://sigproc.sourceforge.net 
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amount of pulse broadening caused by the size of the pre¬ 
vious DM step and then determines the next trial DM. A 
total pulse width smearing due to the DM step in compar¬ 
ison to the value at the last DM is chosen to be 25% (see, 
§ Appendix). The total number of DMs searched was 249, as 
chosen optimally by this program. For the DM steps used 
in our analysis, the average S/N loss is ~ 1.5% for DMs 
< 2000 cm“® pc and for DMs between 2000—10000 cm“® pc, 
the average S/N loss is ^ 2.5%, calculate d using Equations 
12 and 13 of ICordes &i McLaughlinl ll2003ll . The dedisperse 
routine uses Equation 1 to calculate the time delays for each 
test DM and applies to frequency channels and the samples 
from each channel are then averaged to form a dedispersed 
time series. In addition to this, it also performs zero-DM 
subtraction (lEatough et al.l [20091 1 on the time series to re¬ 
move any RFI at zero DM. 


2.2.2 Single-pulse search 

Each dedispersed time series corresponding to a particu¬ 
lar trial DM was searched for transient events of different 
widths via the matched filtering technique for top-hat pulses 
implemented in the program seek. This simple algorithm, 
which is an implementation of the method described in 
ICordes fc McLaughlinl il2003l ~l. saves individual events that 
deviate by five standard deviations from the mean of the 
time series. A number of adjacent samples are added to 
search for pulses of different widths. Each time series was 
smoothed 15 times, corresponding to a maximum smooth¬ 
ing of 2^® times the sampling interval, i.e. pulse widths 
out to 4.096 s. If a pulse is detected in more than one of 
the smoothed times series, only the hig h est S /N value is 
recorded. Following iBurke-Spolaor et al.l (l201lll . the num¬ 
ber of independent trials 


N = Nnm ^ 

3=0 




2J 


( 2 ) 


where Adm is the number of DM trials, Asamp is the number 
of samples in each dedispersed time series, and /max corre¬ 
sponds to number of matched-filter widths used. We find 
A = 1.1 X 10® for our observation. For m o re det ails of this 
sear ch method, seelCordes k, McLaughlinl (l2003ll . 

iKeane fc Petrofll ( 2015l l have assessed the performance 
of the search algorithms commonly being used to dis¬ 
cover FRBs. The authors point out some important con¬ 
cerns where sensitivity to FRBs is often unnecessarily re¬ 
duced and that the single-pulse search routines within 
the dedisperse_all and seek packages were less efficient 
compared to Heim dall and destroy packages (see, e.g., 
IPetroff et al.ll2()14l l. Based on the si mulations which used an 
injected signal of known strength, iKeane fc Petrofll (120151 1 
demonstrated that the recovered S/N for dedisperse_all 
and seek was a function of pulse phase. We apply these 
results into our analysis, as described in § 4. 


2.2.3 Detection criteria 

These dispersed pulses can be displayed graphically using 
the plotpulses program and the resulting single-pulse plot 
can be seen in the left panel of Fig. [T] In the DM histogram, 
we are able to detect sources that emit many weak pulses. 


Such sources might not be detectable in the bottom plot. 
Conversely, the sources that emit only a few strong pulses 
may only be detectable in the bottom plot. Often a peak in 
the DM histogram at low DMs is seen which is indicative of 
RFI. In some cases, RFI is seen at all DMs for a certain time 
range which can further limit our ability to detect a tran¬ 
sient event. This would happen if there is a strong source of 
RFI that causes the receiver to saturate or some other local 
interference that is so strong that it shows up at all DMs for 
that time range. The single-pulse plots thus obtained were 
inspected manually (e.g. by searching for a well-defined peak 
in the S/N versus DM plot, see Fig. [T]). For manual inspec¬ 
tion, we restricted the S/N threshold to six in order to keep 
the number of potential can didates at a manageab le level. 

The event detected by iLorimer et al.l (l2007l l was de¬ 
tected simultaneously in three beams and all of the other 
FRBs were detected in only a single beam. A strong signal 
appearing in all 13 beams simultaneously showing a disper¬ 
sive delay is considered to have a terrestrial or ins trumental 
origin (called pervtons. iBurke-Soolaor et ahllioilll . The ori¬ 
gin of these events has been recently identified as coming 
from a microwave oven when its door is opened prematurely 
and if the telescope is at an appropriate relative angle at 
that time (IPetroff et al.ll2015blj . We did not detect any such 
events in our analysis. We found a number of bursts with 
high S/N in all 13 beams but no dispersive delay in the fre¬ 
quency versus time plot. These events are sources of RFI 
which have near earth origin and are only active for a brief 
period of time. Nearly 10% of the data show S/N greater 
than six in more than five beams but do not show dispersive 
delay. We did not consider these candidates further in our 
analysis. 


2 . 2.4 Manual inspection 

The diagnostic plots obtained after applying all the detec¬ 
tion criteria were manually examined to look for a strong 
signal corresponding to a peak in the S/N versus DM plot. 
Such candidates are shortlisted and the detection is con- 
hrmed if the signal shows a sweep from high to low frequency 
across the observing bandwidth in the frequency versus time 
plot. Brightness was not a criterion for being shortlisted and 
some candidates were confirmed as pulsars despite not be¬ 
ing detectable in this plot. For FRBs, if a dispersive sweep 
is absent then it is deemed to be a “false detection” (see, 
e.g.. Fig. [5]). A burst-like event was seen in one of the beams 
(see a bright pulse in the lower plot of Fig. 0 and a corre¬ 
sponding peak in the S/N versus DM plot in the upper right 
plot) but no dispersed signal was detected in the frequency 
versus time plot (see right panel of Fig. [21). 


3 RESULTS 

As summarized in Table[T] our single-pulse search resulted in 
the detection of 20 of the 42 pulsars detec ted in the original 
periodicity search bv iBurg av et al. (I 2 OO 6 II. I n addition, one 
RRAT not reported by Burgav et al. ( 2006l l was detected. 
Th e discovery of this source, P SR J0410—3113, was reported 
by iBurke-Spolaor et al.l ([SoTJ) during a single-pulse search 
on the data obtained from the high-latitude HTRU survey. 
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4 A. Rane et al. 


Table 1. All pulsars detected and discovered in the PH survey. Columns 1 to 5 report t he pulsar name, Galactic longitude and latitude, 
spin period, and DM, all obtained from the ATNF pulsar catalog llManchester et alj2005lV For those that were detected in our single-pulse 
search method, columns 6, 7, 8, and 9 report the DM obtained in this ana lysis, peak S/N from single-pulse search, number of pulses, 
and width. Column 10 lists S/N f rom the periodicity s e arch obtained from iBurgav et alJ ll200d) and column 11 lists the intermittency 
measure.The RRAT discovered bv iBurke-Spolaor et al. Hm3) and re-detected in our analysis is denoted by *. 


Name 

PSR 

1 

n 

b 

n 

P 

(ms) 

DM 

(pc cm“^) 

DMohs 
(pc cm“^) 

{S/N)sp 

-^pulses 

IF 

(ms) 

(S/Af)per 

r 

J0343-3000 

227.76 

-52.34 

2597.02 

20.2 

22.3 

21.5 

22 

2.0 

42.7 

0.50 

J0410-3113* 

253.47 

-41.95 

1837.00 

9.9 

9.9 

9.2 

1 

4.0 

- 

- 

J0437-4715 

253.47 

-41.95 

5.70 

2.6 

2.6 

12.8 

4303 

0.3 

603.3 

0.02 

J0448-2749 

228.43 

-37.91 

450.40 

26.2 

26.4 

10.2 

6 

4.0 

28.5 

0.38 

J0520-2553 

228.51 

-30.53 

241.60 

33.7 

- 

- 

- 

- 

34.0 

- 

J0610-2100 

227.75 

-18.18 

3.86 

60.6 

_ 

_ 

_ 

— 

10.1 

— 

J0630-2834 

237.03 

-16.75 

1244.40 

34.5 

35.3 

20.2 

63 

16.0 

277.4 

0.07 

J0633-2015 

229.33 

-12.95 

3253.21 

90.7 

89.7 

11.0 

3 

8.0 

16.4 

0.67 

J0636-4549 

254.55 

-21.55 

1984.59 

26.3 

- 

- 

- 

- 

11.2 

- 

J0656-2228 

233.66 

-8.98 

1224.75 

32.4 

31.1 

6.9 

29 

0.3 

21.0 

0.33 

J0719-2545 

238.93 

-5.83 

974.72 

253.9 

_ 

_ 

_ 

— 

20.7 

— 

J0726-2612 

240.08 

-4.64 

3442.31 

69.4 

68.8 

16.8 

9 

4.0 

15.1 

1.11 

J0729-1448 

230.46 

1.44 

251.60 

92.3 

- 

- 

- 

- 

32.5 

- 

J0729-1836 

233.83 

-0.33 

510.10 

61.2 

61.0 

7.4 

26 

0.3 

59.3 

0.13 

J0737-3039A 

245.24 

-4.50 

22.70 

48.9 

- 

- 

- 

- 

18.7 

- 

J0737-3039B 

245.24 

-4.50 

2773.46 

48.9 

_ 

_ 

_ 

— 

_ 

— 

J0738-4042 

254.27 

-9.1 

374.90 

160.8 

161.3 

20.8 

644 

2.0 

542.1 

0.04 

J0742-2822 

243.85 

-2.43 

166.70 

73.8 

73.0 

18.3 

2130 

2.0 

63.1 

0.29 

J0746-4529 

259.20 

-10.10 

2791.03 

134.6 

- 

- 

- 

- 

10.3 

- 

J0749-4247 

257.14 

-8.33 

1095.40 

104.5 

- 

- 

- 

- 

17.6 

- 

J0758-1528 

234.54 

7.24 

682.20 

63.3 

62.9 

13.4 

67 

8.0 

99.0 

0.14 

J0818-3232 

251.36 

1.87 

2161.26 

131.8 

- 

- 

- 

- 

27.3 

- 

J0820-1350 

235.96 

12.61 

1238.10 

40.9 

39.9 

19.0 

28 

16.0 

112.5 

0.17 

J0820-3921 

257.26 

-1.58 

1073.57 

179.4 

- 

- 

- 

- 

13.4 

- 

J0820-4114 

258.82 

-2.72 

545.40 

113.4 

- 

- 

- 

- 

42.6 

- 

J0821-4221 

259.83 

-3.14 

396.73 

270.6 

_ 

_ 

_ 

— 

10.8 

— 

J0823+0159 

222.06 

21.26 

864.80 

23.7 

22.3 

15.8 

54 

8.0 

13.1 

1.21 

J0828-3417 

254.04 

2.58 

1848.90 

52.2 

52.5 

13.3 

5 


49.2 

0.27 

J0835-3707 

257.15 

2.00 

541.40 

112.3 

113.5 

6.9 

1 

8.0 

33.2 

0.21 

J0837-4135 

260.98 

-0.32 

751.60 

147.3 

147.8 

23.2 

303 

4.0 

152.2 

0.15 

J0838-2621 

248.81 

8.98 

308.58 

116.9 

_ 

_ 

_ 

— 

9.6 

— 

J0843+0719 

219.40 

28.22 

1365.86 

36.6 

- 

- 

- 

- 

15.9 

- 

J0846-3533 

257.26 

4.72 

1116.00 

94.1 

- 

- 

- 

- 

103.1 

- 

J0855-3331 

256.92 

7.53 

1267.50 

86.6 

89.7 

9.9 

12 

8.0 

42.5 

0.23 

J0900-3144 

256.16 

9.49 

11.11 

75.7 

- 

- 

- 

- 

20.8 

- 

J0908-1739 

246.19 

19.86 

401.60 

15.8 

_ 

_ 

_ 

— 

23.9 

— 

J0922+0638 

225.48 

36.40 

430.60 

27.3 

26.4 

18.2 

195 

4.0 

96.8 

0.19 

J0944-1354 

249.20 

28.86 

570.20 

12.4 

- 

- 

- 

- 

43.1 

- 

J0953+0755 

228.97 

43.71 

253.00 

2.9 

3.0 

14.0 

657 

2.0 

128.7 

0.11 

J1022+1001 

231.86 

51.11 

16.40 

10.3 

10.4 

9.6 

266 

0.5 

318.7 

0.03 

J1024-0719 

251.77 

40.53 

5.10 

6.4 

_ 

_ 

_ 

— 

33.4 

— 

J1038+0032 

247.15 

48.47 

28.85 

26.6 

- 

- 

- 

- 

18.1 

- 

J1046+0304 

246.48 

51.71 

326.20 

25.3 

28.2 

9.8 

1 

8.0 

25.5 

0.38 
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Figure 1. Example output of one single-pulse search processing pipeline, showing bright single pulses from PSR J0837—4135 around a 
DM value of 147 cm“® pc. The top left panel shows the S/N distribution of the detected pulses, number of pulses versus trial DM (in 
top center), and S/N as a function of trial DM. The lower plot shows S/N of events as a function of time and trial DM. The size of the 
circles is linearly proportional to the S/N of each pulse. The dispersive delay in the frequency versus time plot is seen on the right for 
250 milliseconds of data. 



Figure 2. The left panel shows a non-astrophysical dispersed burst as seen in the lower plot and a corresponding peak in the S/N versus 
DM plot. The non-dedispersed data is then plotted corresponding to time of the peak for one second as seen in the right panel here, 
confirming that it is a false detection. 


In our analysis, only one pulse was detected in this observa¬ 
tion for PSR J0410—31 lib consistent with the non-detection 
in the periodicity search q 

3.1 Energy measurements 

Integrated profiles for each pulsar were obtained by folding 
the dedispersed time series using ephemerides avail able from 
the ATNF pulsar catalog llManchester et ^l2005l ). To con- 
st ruct the pulse en ergy histograms, the procedure described 
in iRitchingsl (Il976l ) is followed. The position and widths of 

® The ATNF pulsar catalog can be accessed online at 
http: //www.atnf.csiro.au/research/pulsar/psrcat 


the on-pulse and off-pulse windows were determined by vi¬ 
sual inspection of the integrated pulsar profile. Baseline es¬ 
timation was done using the off-pulse bins, and was sub¬ 
tracted from the data for each pulse. The total energy for 
each pulse in the on-pulse window was calculated and is 
scaled to account for different widths of on and off windows.. 
The on-pulse data were taken in blocks of about 100 pulses. 
The on-pulse and off-pulse normalized energies, E/E, were 
calculated for each block by dividing on-pulse and off-pulse 
energies within the block by the mean on-pulse energy of 
that block to account for variations due to interstellar scin¬ 
tillation. The pulse energy histograms were constructed for 
10 of the total 17 re-detected pulsars, as shown in Fig. |3l 
based on the number of pulses detected during the total ob- 
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10 



1 10 100 1000 10000 
Number of pulses 

Figure 4. Intermittency measure for each pulsar detection in our 
survey. 


servation time. We did not create energy distribution plots 
for sources with fewer than 100 pulses per observation. 

For some pulsars, the on-pulse distribution peaks 
at the mean energy, as seen for PSR J0742—2822, 
PSR J0630-2834, PSR J0837-4135, and PSR J0922+0638. 
In the case of PSR J0742—2822, the energy histograms sep¬ 
arate out clearly with no zero energy excess in the ON 
histogram, suggesting that this pulsar does not null. For 
other pulsars (PSR J0448—2749, PSR J0828—3417, and 
PSR J1046+0304), the histograms overlap just because the 
S/N is low and not necessarily because the pulsar actually 
nulls. But it is important to note here that the statistics are 
limited by the relatively short observation time in the sur¬ 
vey. So, we did not fit any Gaussians to the histograms since 
the errors would be large. Even if the pulsar is nulling, there 
are insufficient pulses to form a distribution for estimation 
of the nulling fraction in our analysis. 


3.2 Intermittency measure 


The two search algorithms used for pulsar searching show 
varying levels of efficiency w hich depend upon the proper - 
ties of eac h particular pulsar (iMcLaughlin &: Corded 12001^ . 
Following iDeneva et al.1 (120091 '). we calculated the intermit¬ 
tency ratio 


(g/Ar)sp 

(5/A^)per 


(3) 


for each pulsar from the S/N value of the single-pulse and 
periodicity searches and r is plotted versus the number of 
periods in Tobs (see Fig.lT)). All the pulsars re-detected in our 
analysis except PSR J0410—3113 were detected in a previous 
periodicity search. The pulsars on the upper left have longer 
periods and the pulsars on the lower right of this plot are 
millisecond pulsars. Pulsars with r > 1 (PSRs J0726—2612 
and J0823-I-0159) are more likely to be detected with single¬ 
pulse searches. PSR J0726—2612, with r~ 1.1, has P = 3.4 s 
and DM = 69 pc cm~® and PSR J0437—4715, with r ~ 0.02, 
has P = 5.7 ms and DM = 2.6 pc cm“®. These ratios are 
higher than one would expect in longer surveys, because of 
the dependence of sensitivity on number of pulses. The inte¬ 
gration time for this survey was relatively short (265 s) and 


we detected 50% of pulsars from single-pulse search. This 
can be compared to the PM survey with long integration 
time (2100 s), in which the single-pulse detections were only 
about 30% of the total p ulsar discoveries from periodicity 
search ilKeane et al.ll2010ll . 


4 EVENT RATE OF FRBS 


Although our analysis did not result in any new FRB detec¬ 
tions, it is important to consider the impact of this null result 
on the FRB event rate, 77. To put our results into context, 
we also consider a number of other surveys at Parkes in the 
calculation below. To constrain 77, we apply a Bayesian ap¬ 
proach which uses the FR B detections in each s urvey. Using 
Bayes’ theorem (see, e.g.. lWalI fc ,lenkinEll2003h . the poste¬ 
rior probability density function for 77, given the detection 
of n pulses 


p(77|n) oc p(n|77)p(77), (4) 


where p(n|77) is the likelihood of getting n detections given 
some 77 and p(77) is the prior on the rate of FRBs t aken 
to be uniform. Taking n = 0 from IPetroff et al] (|2014h and 
assuming the counting of these rare events as a Poisson pro¬ 
cess, the likelihood function 


P(0|77) 


{TZApTpf exp{-fRTpAp) 
0 ! 


(5) 


where Ap = 4449 deg^ and Tp = 540 sec is the total 
observation time for each pointing. The subscript P cor¬ 
responds to the values for Parkes. The numerical factor 
/ = (86400 X 41253)“^ is inserted to compute 77 in units 
of bursts day“^ sky“^. The posterior probability density of 
the FRB event rate can be computed as 


p(77|0) = i7i77° exp{-fRTpAp), (6) 


where Ki is a normalizing constant which ensures that the 
above expression integrates over all values of 77 to unity. 
Integrating this function, we find that the mean FRB rate 
based on zero FRBs in the HTRU mid-latitude survey is 
77 = 0.22j/g'2i X 10® FRBs day“® sky“®, where the uncer¬ 
tainties represent the 99% confidence interval. To include the 
FRB detections and null results from subsequent surveys, a 
similar calculation is performed to determine the FRB event 
rates for individual surveys as listed in Table 2. The surveys 
were processed using different search algorithms and there- 
fore the rates need to be modified following the results of 
iKeane fc Petrofil (l2015l 'l in order to combine the individual 
ev ent rates. As per the r esponse curve shown in Figure 1 
of iKeane fc Petro3 (l2015l ~). we get the corrected S/N from 
the use of these algorithms by defining the efficiency factors, 
?7ddaii and r^seek- For dedisperse_all, we find 

S/Adedisperse_all ^ 

Jiddaii — e/j\7 — 0.83. (7) 

*5/-^'^max 


and, for seek, 


_ S/Nseek 

“ SjN^ 


0.89. 


( 8 ) 


The FRB searches we consider in this analysis are the re¬ 
sults of_jn 2 roc^sin£_ffie^mkes Multibeam (PMPS) sur¬ 
vey dKeane et al.l (l2010l . [^Olll l. used destroy), the PH sur- 
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vey (this paper, used seek), the S winburne intermediat e 
and high-lati t ude (SWIN) surveys (lEdwards et alJ (2001), 
Ijacobv et alJ (120091 '). iBurke-Spolaor fc BannistCT 1 2014 ). 
used dedisp erse_all), the HTRU hiRh-latitude (H TRU 
high) survey (iThornton et al.l (l2013l ). [ThorntorJ ll201i^ . used 
dedisp erse_all), and the HTRU mid-latitude (HTRU mid) 
survey (IPetroff et al.l ll2015ah . used Heimdall). In addition 
to this, the HTRU surveys were carried out using the digital 
back-end, Berkeley-Princeton-Swinburne Recorder (BPSR), 
whereas the older Parkes surveys used the analogue filter- 
bank (AFB). In order to compare all Parkes surveys, the 
digitization loss factors (/3) depending on the back-end used 
need to be considered in our analysis. For AFB, 0 — 1.25 
and for BPSR, /3 = 1.07 dKouwenhoven fc Voutel[200l]) . We 
then insert the efficiency factors into the minimum fluence 
equation calculated for each processed survey 


Fmin — 


13 T,y,iS/N)„ 


G risott ^NpAf 


Vw, 


(9) 


where Tays is the system temperature, G = 0.66 K/Jy is 
the telescope gain, Np is number of polarizations. A/ is the 
bandwidth, riaoit is either riaeek or r/ddaii depending the sur¬ 
vey, and w is the pulse width equal to the sampling time 
when calculating the minimum fluence. The HTRU mid¬ 
latitude survey has the lowest minimum fluence as can be 
seen in Table [2l The event rates of other surveys are scaled 
to this lowest which is the simplest model assum¬ 

ing a uniform distribution of standard-candle FRBs in Eu¬ 
clidean geometry. The modified event rates are shown in 
Fig. O Combining all these individual rate estimates, Fig. [6] 
shows our current best estimate of TZ. The mean FRB rate 
from this distribution is d.dl®’j x 10® FRBs day“® sky“® for 
sources with a fluence above 4.0 Jy ms at 1.4 GHz, where 
the uncertainties represent a 99% confidence interval. To 
demonstrate that this rate is consistent with all the surveys 
considered here, we list in Table [2] the predicted upper and 
lower bounds on the number of FRBs expected in each sur¬ 
vey which use these 99% confidence intervals on TZ and scale 
it back to each survey’s fluence limit. In addition, we make 
predictions for future analyses of the Perseus Arm (PA) pul¬ 
sar s urvey ( Burgav et al.ll201lj ) and HTRU low-latitude sur¬ 
vey dThornton 20131 ). 

The event rate is estimated assuming FRBs are uni¬ 
formly distributed on the sky. We now consider the impact 
of this assumption. The Galactic effects such as dispersion 
in the ISM, scattering in the ISM, scintillation, and sky tem¬ 
perature can limit the sens i tivity of a survey. Following the 
discussion in IPetroff et al.l (l2014l ) about decreased sensitiv¬ 
ity to FRBs at |b| ^ 15°, we compare the sensitivity based 
on sky coverage for the Parkes surveys considered in our 
calculations. The non-Galactic DM contribution at high lat¬ 
itudes range between 520 and 1070 cm~® pc. The Galac¬ 
tic dispersion at these latitudes is only about 50 cm“® pc, 
whereas the average Galactic dispersion at intermediate lat¬ 
itudes and near the Galactic center are 380 cm“® pc and 
1780 cm“® pc respectively. The FRB pulses with an ad¬ 
ditional non-Galactic DM contribution at these lower lat¬ 
itudes would still be recovered in the surveys considered 
above as they have been searched to a sufficient DM, with 
the maximum trial DM in each of these surveys being 
> 2000 cm“® pc. 

The average sky temperature values for four of these 
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Figure 5. The modified FRB event rates based on individual 
survey results and corrected for detection sensitivity based on 
the search algorithms and the backend used. 



Surveys 


surveys range betwe en 0.85 K to 3.18 K and about 6.14 K for 
the PMPS survey dBurke-Spolaor fc Bannist^ l2014l) . Sky 
temperature is therefore not a significant factor in limiting 
the sensitivity when comparing these surveys. 

FRBs discovered so far (except FRB 110220) show 
few effects of scat tering (see, e.g., lUorimer et al. IH). 
IPetroff et al.l ll2014l ) determined that more than 85% of sur¬ 
vey pointings in the intermediate latitude survey are still 
sensitive to FRB signals even in the presence of strong scat¬ 
tering in the ISM near the Galactic center. This percentage 
will differ for AFB surveys as the number of channels, sam¬ 
pling time, and digitization factors are different but is still 
small compared to overall uncertainty in the FRB event rate 
at t his point. _ 

IPetroff et al.l (l2014l) analyzed how the combination of 
these effects limits sensitivity for survey pointings and de¬ 
termined that a simu lated FRB pulse with properties sim¬ 
ilar to the FRBs in iThorntonI (120131 ) falls below the de¬ 
tection threshold in only 14% of all intermediate latitude 
pointings assuming no scattering in the ISM. We combine 
the individual rates since the sensitivity variations between 
AFB and BPSR ar e within the uncertainty of rate of FRBs. 
IPetroff et al.l (l2014h argue that since the percentage of point¬ 
ings no longer sensitive to FRB pulses decreases to 14%, the 
null result is discrepant with the original predictions based 
on a higher event rate. Although, if the event rate of FRBs is 
much lower than previous estimates, it still explains the null 
result at intermediate latitudes (Table 2). Therefore, we ar¬ 
gue that the lack of detections at intermediate latitudes and 
the null result in our analysis of a high-latitude survey are 
likely to be due to a lower FRB rate and does not necessarily 
imply a dearth of FRBs at intermediate latitudes. 


5 CONCLUSIONS 

We have presented the results of a single-pulse search of the 
PH pulsar survey. W e re-detected 20 of th e previously known 
pulsars reported by iBurgav et al.l (l2006l) and detected one 
RRAT. Out of these 21 re-detections, we have constructed 
pulse energy histograms for 10 pulsars for which the ob¬ 
served number of pulses was more than 100 . For PSR 
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Table 2. Survey parameters, FRB detections, event rates and predicted numbers of FRBs for multibeam pulsar surveys carried out at 
Parkes. The survey abbreviations are listed in the text. For each survey, we provide the dwell time (T), area covered (A), bandwidth 
(A/), digitization factor (/3), efficiency factor [r]), sampling time (taamp), system temperature (Taya), minimum fiuence (Tmin), number 
of FRBs found (W^prbs)! estimated FRB event rate {TZ), and the modified FRB event rate (T^mod)- Using the combined FRB event rate, 
we give the number of FRBs expected to be found (A^expected) i 3-® 3S lower and upper bounds on this expected number (A^io^ver 3nd 

upper). 


Parameter 

HTRU High 

HTRU Mid 

HTRU Low 

SWIN 

PA 

PMPS 

PH 

Dwell time T (s) 

270 

540 

4300 

265 

2100 

2100 

265 

Area A (deg^) 

21671.6 

4448.8 

769.5 

7264.3 

599.3 

2097.4 

3979.5 

Bandwidth A/ (MHz) 

400 

400 

400 

288 

288 

288 

288 

Digitization /3 

1.07 

1.07 

1.07 

1.25 

1.25 

1.25 

1.25 

Efficiency r^goft 

0.83 

1.0 

1.0 

0.83 

1.0 

1.0 

0.75 

Sampling time tsamp (ms) 

0.064 

0.064 

0.064 

0.125 

0.125 

0.25 

0.125 

System temperature Tgys (K) 

29.0 

30.5 

35.5 

29.5 

29.2 

33.6 

29.0 

Fiuence Fmin (Jy ms) 

4.6 

4.0 

4.6 

8.9 

7.3 

11.9 

8.2 

NfRBs 

5 

0 

- 

1 

- 

1 

0 

TZ (10® FRBs d-l aky-l) 

7.9ti« 

0 22^^'^ 
'-'•^^-0.21 

- 

1 

-*-•^-1.8 

- 

0.811^;® 

^•^-0.4 

TZ^od (10® FRBs d~® sky“®) 

9.7li®2° 


- 

’^•^-6.2 

- 


1 c+33.7 
^•^-1.4 

-^expected 

2 

3 

3 

1 

1 

1 

0 

-^upper 

5 

6 

7 

2 

1 

2 

1 

-^lower 

1 

1 

1 

0 

0 

0 

0 



Figure 6. The combined posterior PDF of the FRB event rate 
taking into account all the published Parkes survey results carried 
out to date and correcting for detection sensitivity based on the 
search algorithms used and the backend used. The all-sky event 
rate for sources with a fiuence above 4.0 Jy ms is 7?. = 4.4)j)3’J X 
10® FRBs day“® sky“® . 


J0742—2822, the on-pulse and off-pulse histograms sepa¬ 
rate out clearly, suggesting it does not null over the ob¬ 
serving span. The observation time, being relatively short, 
is insufficient for estimation of the nulling fraction. We 
demonstrated, using a Bayesian approach that the lack of 
FRB detections, and detection rates in other surveys, is 
consistent with an all-sky FRB event rate TZ = x 

10® FRBs day~^ sky”*^, for sources with a fiuence above 
4.0 Jy ms at 1.4 GHz, where the uncertainties represent 
a 99% confidence interval. This event rate takes into ac¬ 
count the decrease in detection sensitivity as a result of 
the search algorithms used in the analysis and the differ¬ 


ent backends used in these surveys. We argue that previ¬ 
ous suggestions of a dearth of FRBs at intermediate Galac¬ 
tic latitudes is likely to be a result of the assumption of 
a higher event rate. Th e revised event rate fou nd here is 
in good agreement with iKeane fc Petroffl (I201!tI ) who esti¬ 
mated TZ ~ 2500 FRBs day“® sky“®. We concur with these 
authors that a larger sample of FRBs is required to more 
meaningfully constrain the rate. 
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7 APPENDIX 

iLevinI (120121 ') describes the the algorithm which accounts for 
the amount of pulse broadening caused by the size of the 
previous DM step and then determines the next trial DM. 
The appropriate step size is computed from the effective 
pulse width (wefi), which is a quadrature sum of the intrinsic 
pulse width (tin), the smearing due to scattering (tscatt), 
dispersion smearing within each frequency channel (Idm), 
sampling time of a particular survey (tsamp), and smearing 
across all frequency channels due to the dispersion measure 
step size (Iadm) 


— \ + Ccatt + ^samp + H 


( 10 ) 


where, for frequencies in GHz and DM in units of cm ^ pc, 


toM = 8300s 


(^) 


( 11 ) 


and 


tADM = 8300s 


/ richanA/ ADM\ 

V 4/3 )■ 


( 12 ) 


Here richan is the number of frequency channels and ADM = 
DM —DM'. Since 4 samples are packed per 64-bit word, the 
equation is divided by a factor of 4. Since the total pulse 
width smearing of 25% is chosen, the pulse broadening frac¬ 
tion due to the DM step is e = 1.25. Then the effective pulse 
width at the new trial DM, with respect to the last trial DM' 
is 


(WeftlDM) — e (w^eff|DM') (13) 

Solving for DM by equating Wes at the new trial value, DM, 
with respect to the last trial value, DM', 


b= DM' + V-a2 b^DM'^ -f a^c -f b^c 
^ a2 -b b2 

(14) 

where 


a - 8300 s , 

r 

(15) 

k k QQnn J^chanA/ ADM 

b — tADM ~ 8300 s „ , 

4 /3 

(16) 

and 


^ ^ Gin 1 ^scatt > ^samp ^DMJ ^samp ^in* 

(17) 


The DM values are obtained using Equation 14 where DM 
step is finely spaced for small values of DM. For higher DM 
values, pulse smearing is greater than the sampling time 
and hence the temporal resolution is decreased giving DM 
channels that are spaced more coarsely. 


© 2015 RAS, MNRAS OOP.IHflOl 













